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Abstract 


Satellite observations of carbon dioxide (CO2) offer novel and distinctive opportunities for 
improving our quantitative understanding of the carbon cycle. Prospective observations include 
those from space-based lidar such as the Active Sensing of CO2 Emissions over Nights, Days, 
and Seasons (ASCENDS) mission. Here we explore the ability of such a mission to detect 
regional changes in CO2 fluxes. We investigate these using three prototypical case studies, 
namely the thawing of permafrost in the Northern High Latitudes, the shifting of fossil fuel 
emissions from Europe to China, and changes in the source/sink characteristics of the Southern 
Ocean. These three scenarios were used to design signal detection studies to investigate the 
ability to detect the unfolding of these scenarios compared to a baseline scenario. Results 
indicate that the ASCENDS mission could detect the types of signals investigated in this study, 
with the caveat that the study is based on some simplifying assumptions. The permafrost thawing 
flux perturbation is readily detectable at a high level of significance. The fossil fuel emission 
detectability is directly related to the strength of the signal and the level of measurement noise. 
For a nominal (lower) fossil fuel emission signal, only the idealized noise-free instrument test 
case produces a clearly detectable signal, while experiments with more realistic noise levels 
capture the signal only in the higher (exaggerated) signal case. For the Southern Ocean scenario, 
differences due to the natural variability in the ENSO climatic mode are primarily detectable as a 
zonal increase. 
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1 Introduction 


Satellite observations of carbon dioxide (CO 2 ) offer novel and distinctive opportunities for 
improving our quantitative understanding of the carbon cycle, which is an important scientific 
and societal challenge with anthropogenic CO 2 emissions still on the rise. Prospective new 
observations include those from space-based lidar such as the Active Sensing of CO 2 Emissions 
over Nights, Days, and Seasons (ASCENDS) mission, which is proposed in "Earth Science and 
Applications from Space: National Imperatives for the Next Decade" [National Research 
Council, 2007] (henceforth referred to as the decadal survey). Notable features of this mission 
include its ability to sample at night and at high latitudes. These conditions are prohibitive to 
passive missions, such as the Greenhouse gases Observing SATellite (GOSAT) [e.g. Kuze et al., 
2009; Yokota et al., 2009] and the Orbiting Carbon Observatory-2 (OCO-2) missions [e.g. Crisp 
et al. 2004], due to their reliance on reflected sunlight. The lidar measurement technique 
proposed for the ASCENDS mission further enables observing through some clouds and aerosols 
[Ehret et al., 2008], which also represent impediments and potential sources of bias for passive 
missions [e.g., Mao and Kawa, 2004]. Extensive instrument design research and development is 
ongoing, and proof of concept and validation studies indicate that ASCENDS will be able to 
provide high-precision, unbiased observations with improved spatial coverage compared to 
passive missions [e.g. Spiers et al, 201 1; Abshire et al., 2010; Kawa et al., 2010]. 

The primary goals of the ASCENDS mission are to address open questions in carbon cycle 
science that focus on the identification of changing source/sink characteristics that are difficult to 
observe using other current or anticipated observations. These goals were first articulated in the 
decadal survey [National Research Council, 2007] and later refined in an ASCENDS mission 
NASA Science Definition and Planning Workshop [ASCENDS Workshop Steering Committee, 
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2008]. They include detecting changes in the northern high latitude sources and sinks, detecting 
changes in Southern Ocean source/sink characteristics, constraining anthropogenic CO 2 
emissions, and increasing our understanding of biospheric carbon dynamics by differentiating 
photosynthetic and respiration fluxes [ASCENDS Workshop Steering Committee, 2008]. CO 2 
fluxes in the northern high latitudes and in the Southern Ocean may change substantially as 
climate evolves, for example, and it is crucial to detect and attribute such changes quickly as they 
could lead to large increases in atmospheric CO 2 concentrations and subsequent shifts in climate 
dynamics [e.g. Canadell etal, 2010]. 

Guided by these stated goals, this study explores the extent to which a space-based lidar 
mission, using the ASCENDS mission concept as a guideline, can indeed contribute to these 
pertinent carbon cycle science questions. We specifically focus on the ability of such a mission 
to detect regional changes in fossil fuel emissions, high latitude CO 2 fluxes, and CO 2 fluxes in 
the Southern Ocean. We investigate these using three prototypical case studies, namely the 
thawing of permafrost in the Northern High Latitudes, the shifting of fossil fuel emissions from 
Europe to the Peoples Republic of China (China), and changes in the source/si nk characteristics 
of the Southern Ocean related to the El Nino Southern Oscillation (ENSO). Realistic flux 
scenarios are defined for each of these prototypical case studies, within the anticipated time 
frame of the ASCENDS mission (i.e., the early to mid-2020s). These flux scenarios, combined 
with a common set of baseline fluxes, form the basis of the presented analysis. 

One can view the experimental setup as a hypothesis testing setup to answer the question if, 
and how, the CO 2 concentration fields resulting from the baseline and perturbation fluxes, as 
observed by an ASCENDS-like mission, are distinguishable. To that end, the three scenarios 
described above are used in Observing System Simulations Experiments (OSSEs) to investigate 
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whether an ASCENDS-like mission will have the ability to identify the changes in atmospheric 
CO 2 distributions associated with the changes in fluxes represented in each scenario. We apply 
the geostatistical mapping approach developed in Hammerling et al. [20 12a, b] to generate global 
CO 2 maps based on ASCENDS-like sampling of the atmospheric CO 2 distribution resulting from 
each flux scenario, and characterize the time required to observe statistically significant changes 
in the inferred global CO 2 distribution, given varying assumptions about measurement 
uncertainty. 

2 Model-Simulated Data 

The study is based on simulated data described below. The study period represents a full 
year of the expected ASCENDS mission data. 

2.1 Baseline CO 2 Atmosphere 

We use the parameterized chemistry and transport model (PCTM) to produce a simulated 
distribution of atmospheric CO 2 variability in space and time \Kawa et al., 2004], based on the 
baseline and perturbation scenarios. Model transport is driven by real-time analyzed 
meteorology from the GEOS-5 Modern-Era Retrospective Analysis for Research and 
Applications (MERRA) [ Rienecker et al., 2011] for 2007. CO 2 surface fluxes for the baseline 
run include terrestrial vegetation physiological processes and biomass burning from CASA- 
GFED3 [Randerson et al., 1996; van der Werfetal., 2010], seasonally-varying climatological 
ocean fluxes from Takahashi et al. [2002], and fossil fuel burning from the CDIAC database 
[Andres et al., 2009]. CASA fluxes are driven by MERRA data, which results in 
meteorologically-driven correspondence in the synoptic variability in the surface fluxes and 
atmospheric transport. The monthly CASA fluxes are downscaled to 3-hourly fluxes in the 
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method of Olsen and Randerson [2004] as described by e.g. Chatterjee et al. [2012] and Shiga et 
al. [2013]. The annual integral for the CASA fluxes is 1.53 Pg/yr representing a net source 
largely due to high respiration in 2007. The average yearly CASA flux for the period 1997 to 
2012 is 0.133 Pg/yr. The CASA fluxes used in this study are available at the North American 
Carbon Data archive (http://nacp-liles.nacarbon.org/nacp-kawa-01/) . PCTM CO 2 output has 
been extensively compared to in situ and remote sensing observations at a wide variety of sites, 
and in most cases the model simulates diurnal to synoptic to seasonal variability with a high 
degree of fidelity [e.g. Law et al., 2008a; Parazoo et al., 2008; Bian et al., 2006; Kawa et al., 
2004]. For the simulations here, the model is run on a 1° x 1.25° latitude/longitude grid with 56 
vertical levels and hourly output. We use 2007 meteorological, cloud and aerosol, and 
reflectivity data for all components employed in the derivation of the simulated ASCENDS CO 2 
observations. 

2.2 Perturbation flux scenarios 

Three case studies are developed based on areas of interest within carbon cycle science that 
are directly relevant to the ASCENDS mission goals, namely the detection of regional changes in 
fossil fuel emissions, high latitude CO 2 fluxes, and changes in CO 2 fluxes in the Southern Ocean. 
They represent quantitatively plausible scenarios of changes in carbon fluxes, henceforth referred 
to as perturbation flux scenarios, that could occur by the early to mid-2020s, the planned launch 
timeframe for ASCENDS. These scenarios are used as prototypical examples of flux patterns 
that give rise to the types of signals the ASCENDS mission endeavors to detect. The 
perturbation fluxes are added to the baseline fluxes within the PCTM modeling framework 
described in section 2. 1 to produce the perturbation CO 2 atmospheres henceforth referred to as 
perturbation runs. 


6 


124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 


2.2.1 Permafrost carbon release 

The permafrost carbon feedback is an amplification of surface warming due to the release of 
CO 2 and methane from thawing permafrost [Zimov et al. 2006]. Permafrost soils in the high 
northern latitudes contain approximately 1700 Gt of carbon in the fonn of frozen organic matter 
[Tarnocai et al. 2009]. Permafrost is perennially frozen ground remaining at or below 0°C for at 
least two consecutive years [Brown et al. 1998] occupying about 24% of the exposed land area in 
the Northern Hemisphere [Zhang et al. 1999]. As temperatures increase in the future and the 
permafrost thaws, the organic material will also thaw and begin to decay, releasing CO 2 and 
methane into the atmosphere. CO 2 and methane emissions from thawing permafrost will amplify 
the warming due to anthropogenic greenhouse gas emissions [Zimov et al. 2006]. 

The permafrost carbon emissions scenario applied here uses projections of CO 2 fluxes from 
thawing permafrost from Schaefer et al. [2011]. Schaefer et al. [2011] use the Simple 
Biosphere/Camegie-Ames-Stanford Approach (SiBCASA) land surface model [Schaefer et al., 
2008] driven by output from several General Circulation Models for the A IB scenario from the 
IPCC Fourth Assessment report [Lemke et al., 2007]. The fluxes are an ensemble mean of 18 
projections from 2002 through 2200. We ran the PCTM model with the extracted fluxes for 2020 
to 2022 and used the 2022 fluxes as perturbation fluxes. The annual integrals for the permafrost 
perturbation fluxes are 0.613 PgC/yr, 0.641 PgC/yr and 0.752 PgC/yr for 2020 to 2022, 
respectively. 

The flux perturbations are concentrated in areas of discontinuous permafrost along the 
southern margins of permafrost regions (Figure 1). In discontinuous permafrost regions, north- 
facing slopes might form pennafrost, while south-facing slopes may not. Permafrost 
temperatures hover just below freezing, making these regions vulnerable to thaw for small 
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increases in atmospheric temperature. Normally, the surface soils in the active layer thaw each 
summer and refreeze each winter. However, as temperatures increase, the thaw depth becomes 
too deep to refreeze in the winter, forming a talik or layer of unfrozen ground above the 
permafrost. The talik allows microbial decay to continue during winter when the surface soils 
are frozen, resulting in year-round fluxes that peak in summer when soil temperatures are highest 
(see Figure A1 in supplementary material). 

The CO 2 concentrations of the baseline run were mean-adjusted to match the annual mean 
of the perturbation run by applying a multiplicative adjustment. This adjustment preserves the 
spatial patterns of the baseline run, while the global difference in concentrations between the 
baseline and perturbation run is zero, so effectively a global flux-neutral scenario. This has been 
done to focus this study on the detectability of changes in spatial patterns rather than detecting 
the mean interannual increase in CO 2 concentrations that results from the strictly positive 
perturbation fluxes over the two years of model spin-up and the investigated year. 

2.2.2 Shift in fossil fuel emissions 

The fossil fuel flux perturbation scenario consists of a shift of fossil fuel emissions from 
Europe to China; a shift that is in directional agreement with recent trends in these regions. 

Fossil fuel emissions from China have increased rapidly over the last decades and China is now 
largest emitter of CO 2 worldwide [ Olivier et al., 2012; Peters et al., 2011]. By comparison, 
fossil fuel emissions from Europe decreased 3% in 201 1 relative to 2010 with an overall decline 
over the last two decades [ Olivier et al., 2012]. We used two magnitudes of emission shift, from 
here on referred to as the “lower” and “higher” signal, representing two points on a continuum of 
possible emission changes around the year 2022. 
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The lower signal represents a 20% decrease of European emissions, with a 12% increase in 
China (Figure 2) that exactly offsets the European decrease. The higher signal includes a 50% 
decrease of emissions in Europe with a corresponding 30% increase in China (Figure 2), and is 
used for illustration purposes only, as a decrease of this size is not expected in Europe within a 
decade. All the percentage changes are in reference to 2007 emission levels, based on the v201 1 
2007 fossil fuel emissions from the CDIAC data base [Andres et al., 2011]. The annual flux 
integrals for the lower and higher signals are 0.228 PgC/yr and 0.571 PgC/yr, respectively. We 
use these two shift settings as examples to draw broader conclusions on the detectability of these 
types of signals as characterized by their the spatial and temporal patterns and their magnitudes. 

The flux perturbations are globally flux neutral, in that European fossil fuel fluxes are 
reduced by a set percentage in each month and the total emissions from China are increased by 
the same mass amount. The decrease and increase is conducted proportionally to the existing 
spatial pattern of the fluxes for each month, thereby preserving the spatial and temporal patterns 
within the European and Chinese fluxes (Figures A2 and A3 in supplementary material). The 
fluxes vary relatively little from month to month, on average +/- 15%. Overall, the signal to be 
detected is a difference in the spatial distribution of CO 2 concentrations, with the global mean 
remaining unchanged. 

2.2.3 Changes in Southern Ocean fluxes 

The Southern Ocean is of special interest to carbon cycle science, because its CO 2 fluxes are 
highly uncertain [ Gruber et al., 2009], it is a region with apparent high sensitivity to climate 
change [Le Quere et al, 2009], and this sensitivity has implications for the region’s future as a 
carbon sink because half of the ocean uptake of anthropogenic CO 2 is estimated to occur there 
[e.g. Le Quere et al., 2009; Meredith et al., 2012]. There is disagreement on the current and 
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future trend of the carbon flux in the Southern Ocean [Le Quere el al., 2009; Law el al., 2008b]. 
The Southern Ocean is also a very sparsely sampled region, where the ability of the ASCENDS 
mission to observe at high latitudes could provide valuable insights. 

Variations in climate modes are a key driver of interannual variability in ocean carbon 
exchange [e.g. Field el al., 2007]. Here we evaluate the extent to which interannual variability 
due to variations in climatic modes, such as the El Nino Southern Oscillation (ENSO), can be 
detected as a reference for addressing potential changes in the si nk /source characteristics of the 
Southern Ocean using satellite observations. In other words, we use ENSO-related variability as 
a prototypical example of the scale of variability to detect. To that end, the years 1977 and 1979 
were chosen as examples of estimated flux patterns, as they represent large differences in ocean 
fluxes due to variations in climatic modes. 

The Southern ocean fluxes for this scenario are based on a hind cast simulation of the 
Community Climate System Model (CCSM) Ocean Biogeochemical Elemental Cycle model as 
described by Doney et al. [2009]. The fluxes were obtained at 1° x 1° spatial and monthly 
temporal resolution. The monthly difference between the Southern ocean flux anomaly for 1977 
and 1979 is used as the perturbation flux for this scenario. Figure 3a and 3c show the average 
flux perturbation for April through June and the full year, respectively. A year-round time series 
of these monthly perturbation fluxes is shown in Figure A4 (supplementary material). The 
magnitude of the perturbation fluxes (1977: 0.186 PgC/yr, 1979: -0.176 PgC/yr) is low relative 
to the other two case studies. In contrast to the other two experiments, the sign of the 
perturbation flux also varies by month and spatially within the region. 
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2.3 Simulated ASCENDS CO 2 observations 


Anticipated ASCENDS sampling and random measurement error characteristics are derived 
from model output and observations in a method similar to that of Kawa et al. [2010] and Kiemle 
at al. [2014]. The Cloud- Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) 
orbital track is used to simulate the expected ASCENDS sampling. Synthetic observations are 
sampled from model output at the nearest time to the satellite over-flight, and interpolated in 
space to the CALIPSO sample locations. A vertical weighting function, appropriate to an 
ASCENDS laser instrument operating at a wavelength near 1.57 um, is applied to the model 
pseudo-data profile to produce column average mixing ratio values [ Ehret et al., 2008]. We 
consider only random errors due to photon counting. Potential bias errors [e.g. Baker et al., 
2010], which could significantly complicate the analysis if correlated with geophysical variables 
of interest (e.g., land/ocean or vegetation cover), are not included. 

CALIPSO measurements of total cloud and aerosol optical depth (OD) are used to calculate 
the ASCENDS laser attenuation. CALIPSO OD data are reported every 5 km (corresponding to 
every 0.7435 s) along track and this forms the basic ASCENDS sample set. Surface lidar 
backscatter ([:!), also needed for error estimation, follows from MODIS measured spectral 
reflectance over land and the glint formulation of Hu et al. [2008] over water using daily 
MERRA 10-m wind speeds. Surface reflectivity over land is interpolated from MODIS (Terra + 
Aqua) 5-km 16-day composite nadir BRDF-adjusted reflectance data (a) at 1.64 pm (band 6), 
which are available every 8 days [Schaafet al., 2002]. Land reflectance is scaled by a factor of 
1.23 to account for the land ‘hot spot’ backscatter effect [Disney et al., 2009], i.e., [3 (sr' 1 ) = 1.23 
a/71. Backscatter values of 0.08 sr" 1 and 0.01 sr’ 1 are used to fill missing areas of MODIS data 
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over land and over snow/ice, respectively, where ice and snow cover is determined from 
MERRA data. 


In order to make our study method applicable for a range of possible CO 2 laser sounder 
instrument implementations, we scale the errors globally to a nominal error value for clear-air 
conditions at Railroad Valley, NV ((] = 0.176, T = 1) and a 10-s (67.2 km) sample integration. 
Thus, a given instrument model can be characterized by its random error at the Railroad Valley 
reference point and the global distribution of errors estimated from OD and p. The individual 
sounding errors at the 5 km CALIPSO resolution are calculated using: 

_ 3.667 * cr ref (1) 

^km ,j3*T 2 *sdf^ 05 ’ 

^ 0.176 ' 

where cr re f is the 10-s reference instrument random error (standard deviation) at Rail Road 
Valley, T the transmittance, self the surface detection frequency, and 0. 176 is the Rail Road 
Valley backscatter reference value at 1.57 pm, which corresponds to one of the potential 
ASCENDS instrument designs [Abshire et al., 2010]. The transmittance is calculated from the 
CALIPSO (OD) using T = e~ OD . 

Soundings with an optical depth greater than 0.3 or where the surface detection frequency 
equals zero are filtered out and considered ’not retrieved'. The surface detection frequency 
equals zero when none of the 1 - km averaged CALIPSO samples in a 5-km average can see a 
ground return, i.e., the clouds/aerosol were too thick to get a return from the ground. For this 
study we a 10-s along-track average as our pseudo-data measurement granule [Kawa et al., 
2010]. Using this setup, the maximum number of soundings constituting one observation is 14. 
The 10-s observation error variances are then calculated by averaging the 5-km error variances 
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within each 10-s time interval and dividing this average by the number of retrieved soundings. 
This setup implies the assumption that retrieval errors between individual soundings are spatially 
and temporally uncorrelated. We include three measurement noise settings in the study. A no- 
measurement noise setting for reference purposes, and medium and high noise settings, which 
use Rail Road Valley 10-s reference instrument random errors (cr ref ) of 0.5 and 1 ppm, 
respectively. Once the measurement error variance has been determined for each location 
following the procedure described above, a random sample from a normal distribution with that 
variance is drawn and added to the PCTM model CO 2 value to define a pseudo-data observation. 
The global mean errors (cr obs ) are 2.1 and 4.2 ppm, respectively, for the medium and high noise 
settings. Figure 4b provides an example of four days of global observations. Different random 
seed numbers are used for the errors in the baseline and the perturbation runs. 

3 Methods 

3.1 Mapping approach 

We use a geostatistical mapping approach [ Hammerling et al., 20 12a, b] to create contiguous 
interpolated maps, i.e. global mapped ("Level 3") products, for the comparison. Satellite CO 2 
observations contain large gaps and high measurement errors such that meaningful spatially- 
comprehensive comparisons at synoptic timescales are often precluded using the observations 
directly. Using gap-filled products makes it possible to conduct synoptic, global comparisons. 
The approach presented in Hammerling et al. [20 12a, b] also yields spatially-explicit 
uncertainties (eq. 7 in Hammerling et al. [2012a]) of the mapped products. These may be lower 
than the uncertainties of the individual observations in areas where the correlation with nearby 
observations can be leveraged, which in turn can facilitate signal detection. 
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The applied mapping methodology yields global mapped CO 2 concentrations with 
uncertainty measures without invoking assumptions about fluxes or atmospheric transport. The 
method leverages the spatial correlation in the atmospheric CO 2 concentration field, 
parameterized as an exponential covariance function with spatially-varying variance and 
correlation range parameters using a moving-window circular domain of 2000 km. These 
parameters are estimated from the observations, and thus not imposed a priori. Methodological 
details are given in Hammerling et al. [2012a, b]. 

Specific to this study, we filter the pseudo-data observations to those with a measurement 
error standard deviation below a certain threshold. This is analogous to imposing a quality 
criterion when delivering remote sensing data instead of making all retrievals available. For the 
medium (high) measurement error setup this threshold is 1.5 (3.0) ppm. This choice represents a 
balance between spatial coverage and robustness of the covariance estimation procedure and was 
determined in a sensitivity analysis (results not shown). For the "no error" setup, which is only 
used as a theoretical best case, we use observations at the same locations as for the medium and 
high measurement setups, but without any noise added. We include this case to isolate any 
potential limitations of the methodology and the spatial coverage from those related to the 
instrument capabilities. 

Mapping CO 2 satellite observations at synoptic time scales makes it possible to capture the 
dynamic behavior of CO 2 in the atmosphere [Hammerling et al., 2012a, b]. Based on 
preliminary studies evaluating mapping performance for different time periods, 4-day periods 
were found to provide the best balance between ensuring good data coverage while also 
capturing synoptic behavior for ASCENDS-like observations. Figure 4 shows an example of a 
4-day (August 1- 4 2007) period. The average modeled CO 2 distribution is shown for reference, 
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and only the observations are used in the subsequent mapping procedure. Each 4-day period is 
mapped independently from other 4-day periods and for the baseline and perturbation runs. For 
January, only six 4-day periods were mapped due to missing CALIPSO data; for all other 
months, seven 4-day periods were mapped for a total of 83 4-day maps for each of the baseline 
and perturbation cases. These mapped fields were then used as input data in the subsequent 
comparison analysis described in Section 3.2 

3.2 Comparison approach 

The detectability of a signal is assessed point-wise for each model grid cell by comparing 
the mapped concentrations from the baseline run and the perturbation run, and determining 
whether the observed differences exceed their associated uncertainties. The uncertainty of the 
difference between two mapped concentration fields, expressed as a variance, is the sum of the 
estimation variances of the baseline, , and the perturbation, a\ , mapped products: 

y base y per 

_2 _ 2 , 2 

G diff Gy, ' Gy 

j j * base J per 

For ease of interpretation and visualization purposes, the comparison results are binned by their 
relative uncertainties, where the absolute value of the difference exceeds one, two or three 
standard deviations, respectively, of the uncertainty of the difference. 

Due to the measurement error added to the observations, together with the sparseness of the 
available observations, individual 4-day maps do not exhibit detectable differences in a statistical 
sense, and the question then becomes to identify a time period over which such 4-day maps must 
be averaged before a significant signal emerges. Under the assumption of temporal 
independence, the uncertainty (expressed as a variance) of the temporal mean is the mean 
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mapping variance of the individual periods divided by the number of periods. The assumption of 
temporal independence was evaluated by conducting temporal variogram analyses for sets of 
mapping errors at randomly selected locations (results not shown), and no compelling indication 
to contradict this assumption was found. 

4 Results and Discussion 

On a high level, one can view the query of detecting atmospheric CO 2 concentration signals 
resulting from flux perturbations as two distinct, if connected, questions. The first question is 
how characteristics of the flux perturbations are translated to, and preserved in, the atmospheric 
CO 2 concentrations, i.e. what is the signature (or signal) of a set of flux perturbations of interest 
in the atmospheric CO 2 concentrations. The second question is how well a given observing 
system, in our case the ASCENDS mission, can capture the presence of this signal. Both of 
these aspects are discussed in the following sections, which are organized by the three 
investigated scenarios. 

4.1 Detectability of permafrost carbon release 

A significant signal can be detected in the case of the anticipated permafrost carbon 
emissions (Figures 5 and A5). The challenge is in capturing longitudinal and latitudinal 
gradients, which can better attribute the increase to the permafrost region, as opposed to just 
detecting a zonal increase. While signal detection is not directly targeted at quantifying carbon 
fluxes, insights on the detectability of spatial gradients are highly relevant for studies targeting 
flux detections, e.g. atmospheric inverse modeling studies. With this in mind, a judicious choice 
of the temporal aggregation periods over which the comparisons are conducted is important. 


16 


341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 


Because of the seasonality of the fluxes in the permafrost carbon release scenario (Figure 
Al), the gradients in the atmospheric CCF distribution are most evident in the months following 
the start of the spring thaw. As a result, averaging over spring/summer months yields a clearer 
identification of the geographic origin of the signal relative to aggregating maps over the full 
year. While the concentration signal is highest around September (Figure Al), or even later in 
the year, when the active layer is deepest, the concentration signals indicative of the spatial 
pattern of the tundra fluxes are most distinct in the late spring/early summer months before the 
effects of atmospheric mixing take over. By August, atmospheric mixing, which occurs rapidly 
in the Arctic, causes the spatial signature of the tundra melting fluxes to be replaced by the 
dominant signal of a zonal increase. Some further evidence of this phenomenon can be observed 
by comparing Figures lb and Id: the 3-month signal retains more of the spatial characteristics, 
whereas the yearly signal represents a zonal increase where the elevated concentrations have 
spread towards the pole. This phenomenon is caused by the specific combination of the temporal 
pattern of the permafrost carbon release and rapid atmospheric mixing in the High Northern 
Latitudes. 

Figures 5 and A5 show a summary of the detection results for the permafrost carbon release 
scenario. While the 3-months results feature comparatively more noise, the recognition of the 
spatial pattern in the significance plots is also improved. Even for the high noise scenario, the 
land origin of the signal is better seen in the 3-month maps relative to the yearly plots. The 
results for the different noise levels are as expected; lower noise provides a more accurate 
mapped concentration field (Figure A5). Overall, the permafrost CO 2 perturbation is detectable 
for both levels of measurement noise considered, and spatial gradients are best detected using 
two to three month aggregation periods in the late spring/early summer. 
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4.2 Detectability of shift in fossil fuel emissions 

The examined fossil fuel emissions perturbations lead to a pronounced spatial signature that 
is localized over Europe and China (Figures 2c, g; 2d, h; A2; A3). This is in contrast to the other 
two experiments (see Figures Id and 3d), where the detailed spatial signatures are largely lost 
and especially the yearly signals represent primarily zonal increases. In addition, the magnitude 
of the lower fossil fuel perturbation signal is very low, which, combined with its small spatial 
extent, renders the weak fossil fuel signal the most difficult to detect among the investigated 
scenarios. This exemplifies the challenge of detecting small and localized flux changes from 
satellite observations. 

The highly localized nature of the signal over Europe and China suggests that changes in 
fluxes would be detectable and attributable to a given region, even if these fluxes were not offset 
by a corresponding shift in emissions within a similar latitudinal range. It is interesting to 
observe the different dispersion patterns of the European and Chinese emissions, especially when 
considering their latitudinal range. The effect on the atmospheric concentrations of the European 
emissions can be observed from equatorial Africa to the Arctic, whereas the range of the Chinese 
emissions is more limited to their originating latitudinal band (see Figure 2h). 

Given the relative lack of seasonality in the fossil fuel perturbation scenarios, averaging 
over longer periods of time leads to better detectability (Figure 6). Although atmospheric 
transport clearly plays a role, the atmospheric signal remains indicative of the source region of 
the perturbation flux throughout the seasons. Figure 6h, for example, shows evidence that both 
the source and downwind regions of the emissions have a significant signature in the 
atmosphere. 
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The effect of varying measurement noise levels on the detectability is again as expected; 
increasing measurement noise leads to decreased significance in the results and requires in turn 
longer averaging periods. For the higher signal, all three noise levels capture the signal in the 
yearly results, which is not the case for the lower signal, where only the no error case clearly 
captures the signal (Figure A6). There is some evidence that a significant signal is detectable in 
the yearly medium and high error measurement noise cases (Figure 6), and given the nature of 
the signal discussed above, the signal is expected to appear more clearly when averaging over 
periods exceeding one year. Overall, these findings imply that ASCENDS can in principle detect 
anthropogenic signal components, but depending on the strength of the signal, detection might 
require multiple years. It is hence feasible that ASCENDS can serve to validate anthropogenic 
emission changes over the course of its mission, but is likely not ideal as the primary monitoring 
tool for such flux changes. 

4.3 Detectability of changes in Southern Ocean fluxes 

The detection of changes in the Southern Ocean source/si nk characteristics is challenging as 
a result of a confluence of different factors. The overall magnitude of the signal in the Southern 
Oceans is rather weak, with the absolute value of the signal never exceeding 0.4 ppm in the 
column. In addition, this scenario involves sub-seasonal and sub-regional-scale flux variability 
that is superimposed on a seasonal pattern in the fluxes (Figure A4). Atmospheric mixing also 
plays a role insofar as it obscures the Southern Ocean as the origin of the signal, as was also 
observed in the permafrost carbon release scenario. However, applying the remedy of using a 
shorter averaging period before atmospheric mixing hides the origin of the signal is not as clear- 
cut for the Southern Ocean scenario, because the overall signal is weaker. Spatial gradients 
associated with the Southern Ocean flux perturbation are most evident in the spring and early 


19 


409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 


summer (Figure 3). Later in the year, although the concentration signal is stronger, the 
concentration increase has spread pole-ward and is less attributable to the Southern Ocean. 

For all measurement noise setups, the yearly results clearly indicate a zonal increase in the 
High Southern Latitudes (Figures 7 and A8). However, it is less clear whether the pattern is 
indicative of the Southern Ocean being the source region within the zonal band. The spatial 
pattern of the 3-month results (Figures 7b) is more indicative of the Southern Ocean as the 
source region, but the significance levels are not very high. The most beneficial approach for the 
Southern Ocean scenario appears to be conducting analyses over periods of multiple lengths, and 
drawing conclusions from the joint picture emerging from these analyses. In summary, 
ASCENDS can detect a Southern Ocean signal representative of differences due to natural 
variability in the ENSO climatic mode. ASCENDS may provide a unique measurement view of 
these regions because the pervasive cloudiness and low sun angles present difficult conditions 
for passive satellite CO 2 sensors. Due to the low magnitude and small-scale variability within 
the fluxes giving rise to the signal, however, attributing the signal to specific ocean regions or 
biogeophysical processes will likely require additional corroborative infonnation. 

5 Conclusions 

This work assesses the degree to which ASCENDS, a planned lidar CO 2 observing satellite 
mission, can contribute to the detection of three types of CO 2 flux change scenarios relevant to 
global carbon cycle science: the release of carbon due to the thawing of permafrost in the 
Northern High Latitudes, the shifting of fossil fuel emissions from Europe to China, and ENSO- 
related changes in the source/si nk characteristics in the Southern Ocean. These three scenarios 
were used to design OSSEs for signal detection studies to investigate if the ASCENDS mission 
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has the ability to detect the unfolding of these scenarios compared to a baseline scenario. Two 
different levels of measurement noise and a no measurement noise reference case were 
investigated. 

This study is based on a number of simplifications. For each scenario, the only flux 
component that is varied is the flux component under investigation, while all other fluxes are 
fixed. In reality, many changes might occur simultaneously and the resulting CO2 concentration 
signal patterns might overlap, which makes signal detection more challenging. We have 
introduced some additional variability by sampling and mapping the baseline concentration field 
rather than assuming a static baseline concentration field in the comparison procedure, however, 
that might not be equivalent to, for example, having misspecified biospheric fluxes. Such 
misspecifications could be aliased with the true signals and misleading signal patterns could 
occur. This could impact the conclusions of this study insofar that it would be more difficult to 
link detectable signals with the underlying change in fluxes. 

The results indicate that the ASCENDS mission can in principle detect the types of signals 
investigated in this study. The permafrost thawing flux perturbation is readily detectable at a 
high level of significance. Spatial gradients, which are of great interest for process attribution, 
were best detected using two or three month aggregation periods in the late spring/early summer. 
For the Southern Ocean scenario, differences due to the natural variability in the ENSO climatic 
mode were primarily detectable as a zonal increase. The relative magnitude of the signal, 
however, is much smaller than the permafrost-thawing signal. Spatial and temporal high- 
frequency changes in the anomaly fluxes produce additional variability in the signal, making 
detection of more detailed gradients than a zonal increase challenging for the Southern Ocean 


21 


453 

454 

455 

456 

457 

458 

459 

460 

461 

462 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 


scenario. Conducting analyses over periods of varying lengths and analyzing them jointly 
provides a possible diagnostic strategy. 

The fossil fuel emission detectability is directly related to the strength of the signal and the 
level of measurement noise. As is true for all scenarios, the effect of varying measurement noise 
levels is as expected: increasing measurement noise levels lead to decreased significance in the 
results and require in turn longer averaging periods. For the nominal (lower) fossil fuel emission 
signal, only the noise-free instrument test produces a clearly detectable signal, while all three 
noise levels capture the higher (exaggerated) signal case. The emergence of a detectable signal 
suggests that averaging over periods longer than the one-year period considered in this study 
would also render signals of the magnitude of the lower fossil fuel emission signal detectable. 

All in all, the expected precision and sampling characteristics of ASCENDS promise to 
substantially enhance our ability to detect variations in CO 2 fluxes and to infonn the mechanisms 
that control them. Future work includes comparing the signal detection performance of 
ASCENDS to passive sensors, which might be employed within the time frame of the 
ASCENDS mission. Additional future work entails a comprehensive study of the effect of 
uncertainties in fluxes other than those defining the signal on the detectability of the signal, for 
example, by using an ensemble of biospheric fluxes to vary the baseline fluxes. 
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Figure 1 : Flux and CO 2 concentration for the permafrost carbon release experiment, (a) 3- 
month average CO 2 flux (“3-month flux”), (b) 3-month average CO 2 concentration (“3-month 
signal”), (c) yearly average CO 2 flux (“Yearly flux”), and (d) yearly average CO 2 concentration 
(“Yearly signal”). The 3 -month period is May through July. The flux is modeled for 2022. The 
negative concentration values in the Southern hemisphere are a result of the global mean 
adjustment. 
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Figure 2: Flux and CO 2 concentration for the fossil fuel experiments. First row: 3-month 
average CO 2 flux (“3-month flux”). Second row: 3-month average CO 2 concentration (“3-month 
signal”). Third row: Yearly average CO 2 flux (“Yearly flux”). Fourth row: Yearly average CO 2 
concentration (“Yearly signal”). The 3-month period is August through September. 
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Figure 3: Flux and CO 2 concentration for the Southern Ocean experiment, (a) 3-month 
average CO 2 flux (“3-month flux”), (b) 3-month average CO 2 concentration (“3-month signal”), 
(c) yearly average CO 2 flux (“Yearly flux”) and (d) yearly average CO 2 concentration (“Yearly 
signal”). The 3-month period is April through June. 
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Figure 4: Mapping results for August 1-4 2007. (a) Modeled CO 2 concentrations (“Model”), 
(b) simulated ASCENDS observations (“Observations”), (c) mapped CO 2 concentrations 
(“Mapped”), and (d) mapping uncertainties (“Uncertainty”) expressed as a standard deviation. 
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Figure 5: Results for the permafrost carbon release experiment for medium measurement 
noise, (a) 3-month mapped CO 2 signal (“3-month mapped”), (b) Significance of the 3-month 
mapped CO 2 signal (“3-month signific.”), (c) Yearly mapped CO 2 signal (“Yearly mapped”), and 
(d) Significance of the yearly mapped CO 2 signal (“Yearly signific.”). The mapped signal is the 
difference between the mapped perturbation CO 2 concentration and the mapped baseline CO 2 
concentration. The significance is the mapped signal divided by the uncertainty of the mapped 
signal. The values are discretized for improved visualization. Yellow, orange and dark red 
(light, medium and dark blue) represent areas where the mapped perturbation concentration is 
larger (smaller) than the mapped baseline concentration by more than one, two or three standard 
deviations, respectively, of the uncertainty of the mapped signal. The 3-month period is May 
through July. 
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Figure 6: Results for the fossil fuel experiments for medium measurement noise. First row: 
3-month mapped CCF signal (“3-month mapped”). Second row: Significance of the 3-month 
mapped CCF signal (“3-month signific.”). Third row: Yearly mapped CCF signal (“Yearly 
mapped”). Fourth row: Significance of the yearly mapped CCF signal (“Yearly signific.”). The 
mapped signal is the difference between the mapped perturbation CCF concentration and the 
mapped baseline CO 2 concentration. The significance is the mapped signal divided by the 
uncertainty of the mapped signal. The values are discretized for improved visualization. 
Yellow, orange and dark red (light, medium and dark blue) represent areas where the mapped 
perturbation concentration is larger (smaller) than the mapped baseline concentration by more 
than one, two or three standard deviations, respectively, of the uncertainty of the mapped signal. 
The 3-month period is August through September. 
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Figure 7: Results for the Southern Ocean experiment for medium measurement noise, (a) 3- 
month mapped CO 2 signal (“3-month mapped”), (b) Significance of the 3-month mapped CO 2 
signal (“3-month signific.”), (c) Yearly mapped CO 2 signal (“Yearly mapped”), and (d) 
Significance of the yearly mapped CO 2 signal (“Yearly signific.”). The mapped signal is the 
difference between the mapped perturbation CO 2 concentration and the mapped baseline CO 2 
concentration. The significance is the mapped signal divided by the uncertainty of the mapped 
signal. The values are discretized for improved visualization. Yellow, orange and dark red 
(light, medium and dark blue) represent areas where the mapped perturbation concentration is 
larger (smaller) than the mapped baseline concentration by more than one, two or three standard 
deviations, respectively, of the uncertainty of the mapped signal. The 3-month period is April 
through June. 
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